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Abstract — A  standard  assumption  in  most  tracking  algorithms, 
like  the  Probabilistic  Data  Association  (PDA)  filter.  Multiple  Hy¬ 
pothesis  Tracker  (MHT)  or  the  Multiframe  Assignment  Tracker 
(MFA),  is  that  a  target  is  detected  at  most  once  in  a  frame 
of  data  used  for  association.  This  one-to-one  assumption  is 
essential  for  correct  measurement-to-track  associations.  When 
this  assumption  is  violated,  the  above  algorithms  treat  the  extra 
detections  as  random  clutter.  When  multiple  detections  from 
the  same  target  fall  within  the  association  gate,  the  PDA  filter 
tries  to  apportion  the  association  probabilities,  but  with  the 
fundamental  assumption  only  one  of  them  is  correct.  The  MFA 
and  the  MHT  algorithms  try  to  spawn  multiple  tracks  to  handle 
the  additional  measurements  from  the  same  target,  assuming  at 
most  one  measurement  came  from  each  target.  Both  of  these 
approaches  have  undesirable  side  effects  since  they  ignore  the 
possibility  of  multiple  detections  from  the  same  target  in  a  scan 
of  data.  Such  multiple  detection  situations  occur  in  multistatic 
tracking  problems. 

In  this  paper,  we  proposed  a  new  Multiple  Detection  Proba¬ 
bilistic  Data  Association  (MD-PDA)  filter  for  tracking  a  target 
when  more  than  one  target  originated  measurement  may  exist 
within  the  validation  gate.  In  the  proposed  MD-PDA,  combina¬ 
torial  association  events  are  formed  to  handle  the  possibility  of 
multiple  measurements  from  the  same  target.  Modified  associa¬ 
tion  probabilities  are  calculated  with  the  explicit  assumption  of 
multiple  detections.  Simulations  are  presented  to  demonstrate  the 
effectiveness  of  the  algorithm  on  a  single  target  tracking  problem 
in  clutter.  Extensions  to  handle  multiple  targets  using  the  Joint 
PDA,  MHT  and  MFA  approaches  are  under  development. 
Keywords:  target  tracking  in  clutter,  multistatic  tracking, 
data  association,  probabilistic  data  association 

I.  Introduction 

Tracking  a  target  in  clutter  where  it  is  unknown  which  of 
the  received  set  of  measurements  is  originated  from  target 
has  been  one  of  the  most  challenging  issues.  In  the  literature, 
several  techniques  have  been  proposed  for  this  data  association 
problem  to  identify  target  originated  measurement  from  a 
clutter  [2],  [4],  [6],  [12],  [13].  Non-Bayesian,  one-to-one 
matching,  hard  decision  oriented  data  association  solutions  are 
the  Nearest  Neighbor  Filter  (NNF)  and  Strongest  Neighbor  Fil¬ 
ter  (SNF)  [2].  As  their  names  imply,  the  NNF  updates  a  track 
with  the  measurement  closest  to  the  predicted  measurement 
among  the  validated  measurements  while  the  SNF  associates 
the  measurement  with  the  strongest  intensity. 


The  aforementioned  data  association  techniques  perform 
well  in  terms  of  computation  and  estimation  accuracy  in  a 
scenario  where  the  target  return  is  very  strong  and  the  false 
alarm  rate  is  low.  With  degraded  observability  and  dense 
clutter  such  approaches  begin  to  fall  short.  Under  such  con¬ 
ditions,  a  more  practical  approach  to  deal  with  measurement 
origin  uncertainty  is  applying  Bayesian  association  techniques. 
One  of  Bayesian  association  approach  is  the  sub-optimal 
Probabilistic  Data  Association  (PDA)  filter  [3],  [4],  [5].  The 
PDA  estimator  avoids  a  hard  association  decision  by  updating 
a  track  with  a  set  of  measurements  and  their  corresponding 
weights.  In  the  PDA  estimator  the  weight  corresponding  to 
each  validated  measurement  is  calculated  by  assessing  all 
possible  measurement-to-track  combinations.  As  a  result,  the 
weight  assigned  to  a  given  measurement  inside  the  validation 
gate  is  the  probability  that  it  came  from  that  track. 

The  PDA  estimator  is  appropriate  for  single  target  tracking 
problem.  The  Multiple  Target  Tracking  (MTT)  brings  more 
challenge  to  the  data  association  as  it  has  to  be  determined 
which  measurement  belongs  to  which  target  besides  iden¬ 
tifying  target  originated  measurement  from  a  clutter.  Thus, 
the  PDA  has  to  be  extended  to  handling  multiple  targets, 
resulting  in  the  joint  PDA  (JPDA)  algorithm  that  can  handle 
tracking  multiple  targets  by  evaluating  the  joint  probabilities 
among  the  tracks  and  the  measurements  [13].  In  addition, 
the  interacting  multiple  model  is  usually  integrated  to  handle 
target  maneuvers  in  IMM-JPDA  [7].  Furthermore,  multiple 
scan  JPDA  and  multipattern  algorithm  have  been  developed 
as  an  extension  to  the  PDA  [9],  [15]. 

An  alternative  optimal  Bayesian  approach  to  MTT  is  the 
multiple  hypothesis  tracker  (MHT)  [6].  MHT  handles  the 
multitarget  tracking  problem  by  forming  multiple  hypotheses 
and  evaluating  the  likelihood  that  there  is  a  target  in  a 
given  sequence  of  measurements.  Although  it  is  an  optimal 
approach,  within  few  steps  it  will  become  computationally 
infeasible.  Hypothesis  pruning  techniques  can  be  applied  to 
the  MHT  approach  for  practical  problems  at  the  expense  of 
optimality  [10].  Another  approach  to  the  MTT  problem  is 
Multiframe  Assignment  Tracker  (MFA)  [8]  that  models  the 
measurement-to-track  association  as  a  constrained  optimiza- 
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tion  problem. 

The  assumption  in  most  tracking  algorithms,  like  the  afore¬ 
mentioned  PDA,  JIPDA,  MHT  or  the  Multiframe  Assignment 
Tracker  (MFA),  is  that  a  target  is  detected  at  most  once  in 
a  frame  of  data  used  for  association  [6],  [8],  [12],  [14],  This 
one-to-one  assumption  is  essential  for  correct  measurement- 
to-track  associations.  When  this  assumption  is  violated,  for 
example,  when  a  target  is  detected  more  than  once  per  scan, 
the  above  algorithms  treat  the  extra  detections  as  random 
clutter  or  tend  to  spawn  multiple  tracks  for  the  same  target. 

Such  multiple-detection  situations  occur  in  multistatic  track¬ 
ing  problems.  When  multiple  detections  from  the  same  target 
fall  within  the  association  gate,  the  PDA  estimator  as  well  as 
its  multitarget  version,  JIPDA,  try  to  apportion  the  association 
probabilities,  but  with  the  fundamental  assumption  only  one 
of  them  is  correct.  The  MFA  and  the  MHT  algorithms  try  to 
spawn  multiple  tracks  to  handle  the  additional  measurements 
from  the  same  target,  due  to  the  basic  assumption  at  most  one 
measurement  came  from  each  target.  Both  of  these  approaches 
have  undesirable  side  effects  since  they  ignore  the  possibility 
of  multiple  detections  from  the  same  target  in  a  scan  of  data. 
A  mechanism  that  accounts  for  the  possibility  of  multiple 
detections  from  the  same  target  should  be  developed  so  that 
all  useful  information  from  the  received  measurements  about 
the  target  state  is  extracted. 

In  order  to  rectify  the  above  short  coming,  a  new  Multiple 
Detection  Probabilistic  Data  Association  (MD-PDA)  is  pre¬ 
sented  in  this  paper.  In  the  proposed  MD-PDA,  combinatorial 
association  events  are  formed  to  handle  the  possibility  of  mul¬ 
tiple  measurements  from  the  same  target.  Multiple  association 
events  are  formed  by  creating  p  out  of  m  combinations  of 
multiple  measurements  to  track  assignment,  where  p  is  the 
number  of  target  originated  measurements  and  to  is  the  total 
number  of  measurements  in  the  validation  gate.  The  number  of 
target  originated  measurements  can  be  used  as  a  known  prior 
to  determine  the  probability  of  detection  condition  on  p,  which 
is  Pdip-  For  each  association  events,  the  probabilities  will  be 
calculated  and  based  on  the  probabilities  a  measurement  or  set 
of  measurements  will  be  associated  to  a  target. 

The  enhanced  capability  of  the  fusing  all  the  information 
from  target  originated  measurements  in  the  proposed  MD-PDA 
will  manifest  if  the  sensor  data  contains  multiple  detection 
in  a  single  frame.  If  the  target  is  detected  only  once  per 
frame  the  MD-PDA  filter  cannot  perform  better  than  original 
PDA.  Simulation  is  done  by  generating  multiple  detection 
measurements  from  a  single  target  observed  in  a  clutter. 
MD-PDA  performance  is  compared  with  the  original  PDA. 
Performance  evaluation  results  show  the  effectiveness,  with 
respect  to  estimation  accuracy,  of  the  proposed  algorithm  as 
a  result  of  taking  the  possibility  of  multiple  detections  into 
account.  However,  the  algorithm  tends  to  take  more  time  due 
to  increased  number  of  association  events.  Similar  extensions 
to  handle  duplicate  detections  in  the  presence  of  multiple 
targets  using  the  JPDA,  MHT  and  MFA  approaches  are  under 
development. 

The  remainder  of  the  paper  is  organized  as  follows.  Sec¬ 


tion  II  discusses  the  multiple  detection  pattern.  Models  for 
combinatorial  events  in  the  presence  of  multiple  detection  is 
presented  in  this  section.  The  new  MD-PDA  filter  is  presented 
in  Section  III  where  theoretical  developments  of  MD-PDA  are 
discussed.  Simulation  results  is  presented  Section  V,  which  is 
based  on  a  single  target  tracking  in  clutter.  Finally,  conclusion 
is  drawn  in  Section  VI. 

II.  Multiple  Detection  Pattern 

When  multiple  detections  from  the  same  target  fall  within 
the  association  gate,  a  measurement  or  set  of  measurements 
might  be  associated  to  a  target.  Data  association  uncertainty 
with  multiple  detection  can  be  resolved  by  generating  a 
multiple  detection  pattern.  The  multiple  detection  pattern  will 
consider  all  possible  events  for  measurement-to-track  associ¬ 
ation. 

Assume  that  the  targets  state  evolves  according  to  a  dynamic 
equation  driven  by  process  noise 

x{k  +  1)  =  F{k)x{k)  +  w{k)  (1) 

and  the  measurement  equation  by 

z(k)  =  H(k)x(k)  +  v(k)  (2) 

where  x(k)  represents  target  state,  F(k)  is  the  system  transi¬ 
tion  matrix  and  H{k)  is  the  measurement  matrix.  w{k)  and 
v(k)  are  white  and  independent  system  and  measurement  noise 
respectively. 

For  to  number  of  measurement  inside  the  validation  gate 
p  out  of  to  association  events  are  evaluated  while  p  runs 
from  one  to  the  maximum  number  of  target  originated 
measurements.  This  association  event  represent  all  possi¬ 
ble  events  from  single  target  originated  measurement  to 
all  of  the  measurements  are  target  originated.  For  exam¬ 
ple,  as  depicted  in  Figure  1,  there  are  four  measurements 
Z2{k),  z3(k),  Zi(k))  in  the  data  frame.  Out  of  four  mea¬ 
surements,  three  of  them  (zi(fc),  22(A)),  Zz(k))  are  inside  the 
validation  gate.  Combinatorial  association  events  are  created 
only  for  those  measurements  that  fall  inside  the  validation  gate. 
The  maximum  number  of  target  originated  measurement  is 
assumed  to  be  pmax  =  3.  Thus  the  possible  events  are: 

•  none  of  the  measurements  is  target  orignated 

-  p  =  0,  =  1 

•  one  of  the  measurements  is  target  originated 

-  p  =  1,  nv  =  comb(3, 1)  =  1,2,3 

-  3  measurement-to-track  association  events 

-  Z\(k)  or  22 (k)  or  23(A))  is  originated  from  a  target 


21,1(fc)  =  2i(fc) 

(3) 

21.2(A))  =  Z2{k) 

(4) 

21.3(A))  =  z3(k ) 

(5) 

*  two  of  the  measurements  are  target  originated 

-  p  =  2,  nv  =  comb(3, 2)  =  1, 2, 3 

-  3  measurements-to-track  association  events 
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Figure  1.  Validation  Gate 


-  zi(k),  z2(k)  or  zi(k),  z3(k)  or  z2(k),  z3(k)  are  orig- 


inated  form  a  target 

22,1  [k)  = 

zi(k) 

_  z2(k) 

(6) 

Z2,2{k)  = 

zi(k) 

_  23 (k)  _ 

(7) 

22,3  (k)  = 

Z2{k) 

_  z3(k)  _ 

(8) 

•  all  of  the  measurements  are  target  originated 

-  ip  =  3,  =  comb(3, 3)  =  1 

-  1  measurements-to-track  association  event 

-  z\(k) ,  z2(k) ,  z3{k)  are  originated  form  a  target 


23,l(fc) 


zi{k) 

Z2(k) 

zz{k) 


(9) 


Accordingly  the  measurement  equation  (2)  for  the  (p,nv) 
event  becomes 


'  H^k)  ' 

’  vi(k) 

(&)  — 

x{k)  + 

_  Hv(k)  _ 

_  vv(k) 

III.  Multiple  Detections  Probability  Data 
Association  (MD-PDA) 

The  approach  of  the  standard  PDAF  is  to  calculate  the 
association  probabilities  for  each  validated  measurement  that 
falls  in  a  gate  around  the  predicted  measurement  at  the 
current  time  to  the  target  of  interest  [2].  When  two  of  the 
measurements  are  target  originated,  the  algorithm  apportion 
the  total  weight  to  both  of  them,  with  the  assumption  that  only 
one  of  them  is  target  originated.  This  is  not  efficient  approach 
especially  when  there  are  false  alarms  in  the  validation  gate. 
The  MD-PDA  algorithm  calculates  the  probability  that  each 
set  of  measurements,  rather  than  a  single  measurement,  is 


attributable  to  the  target  of  interest.  The  set  of  measurements 
candidate  for  evaluation  are  generated  from  multiple  detection 
pattern  discussed  above.  This  probabilistic  (Bayesian)  infor¬ 
mation  based  on  the  candidate  set  of  measurements  is  used  in 
a  tracking  filter,  that  updates  the  target  states. 

A.  Assumptions 

The  following  assumptions  are  made  for  the  MD-PDA  filter 

•  Among  the  validated  measurements,  a  measurement  or 
set  of  measurements  can  originate  from  a  target. 

•  The  target  detections  occur  independently  over  time  with 
known  probabilities. 

•  Clutter  is  uniform/Possion  distributed  within  the  measure¬ 
ment  validation  gate. 

•  There  is  only  one  target  of  interest  whose  state  evolves 
according  to  a  dynamic  equation  driven  by  process  noise 
as  stated  in  (1). 

•  Track  has  been  initiated. 

At  each  time  k,  the  MD-PDA  algorithm  runs  through  the 
following  steps. 

B.  Gating 

A  validation  gate  is  set  up  for  each  time  step  to  determine 
the  candidate  measurements  for  association.  The  validation 
gate  is  an  ellipse  [2]  given  by 

y(fc, 7)  =  {z  :  [z  —  z(k\k  —  l)],S'(fc)_1 

[z-z(fc|fc-l)]<7>  (11) 

where  7  is  the  gate  threshold  and  S(k)  is  the  innovation 
covariance  corresponding  to  the  measurement  given  by 

S(k)  =  H{k)P{k\k  -  1  )H{k)'  +  R(k)  (12) 

The  volume  is  thus  given  by 

V(k)  =  CnMS(k)  |1/2  (13) 

=  cn^|(S(fc)|1/2  (14) 


where  nz  is  the  dimension  of  the  measurement  and  the  coeffi¬ 
cient  cUz  is  the  volume  of  the  n, -dimensional  unit  hypersphere 

(ci  =  2,  c2  =  7 r,  c3  =  47t/3,  etc.). 


C.  MD-PDA  approach 

The  latest  set  of  validated  measurements  is  denoted  as 

m  =  {*(*)}£?>  (is) 

where  Zi{k)  is  the  ith  validated  measurement  and  m(k)  is  the 
number  of  measurements  in  the  validation  region  at  time  k. 
The  cumulative  set  of  measurements  up  to  time  step  k  is 

Zk  =  (16) 


For  the  association  events 

(ip  out  of  m(k)  are  target  originated  ) 


(k)  — 


tlcp  1)  ***)  ttfpxyi  (^0 


(none  of  the  measurements  is  target  originated) 

(17) 


nv  =  0 
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where  cvm(fc)  is  ip  combinations  out  of  m(k)  measurements 
given  by 

<*»(*>  =  (  ^  )  (18) 

The  number  of  association  events  grow  very  fast  for  ip  >  2. 
The  expected  number  of  target  originated  measurement  can  be 
used  as  a  priori  to  reduce  the  number  of  events.  Applying  the 
total  probability  theorem  with  respect  to  the  above  events,  the 
conditional  mean  of  the  state  at  time  k  is  given  as 

x(k\k)  =  E{x{k)\Zk) 

m{k)  ciprn(k ) 

-EE  E(mK,nv(k),Zk)p(9v,nv(k)\Zk) 

<p= 0  n<p  =  l 

m(k)  clfirn(k) 

=  E  E  ^  (k\k)P<p,nv{k)  (19) 

lp=  0  nv  =  l 

where  (k\k)  is  the  updated  state  which  is  conditioned 

on  the  event  that  the  (tp,  n^)  set  of  measurements  are  correct. 
Here  the  association  probability,  Pv,nv{k),  is  the  conditional 
provability  of  the  event. 

P<p,nv(k)  (xp(ev,nv(k)\Zk)  (20) 

The  estimate  conditioned  on  ntk  combination  of  ip  measure¬ 
ments  being  correct  is 

x<p,nv{k\k)  =  x(k\k  -  1)  +  WVtnv(k)plf>tnv(k) 

where  the  corresponding  innovation  is 

(. z(k )  —  z(k\k  —  1))' 

Vv,nv{k)=  \  (21) 

(. z(k )  —  z(k\k  —  1))' 
and  the  Kalman  gain  W(k )  is 

Wv,n^k)  =  P(k\k  -  1  )H^nv(k),S^nv(k)~1  (22) 

Here 

Sv>nv(k)  =  Hv,n(k)P(k\k-  l)Hv,nv(k)' +  Ry^ik)  (23) 

'  H{k)  - 

Hv,nv{k)  =  :  (24) 

_  H(k)  _ 

'  R{k)  0  ...  0 

0  R(k)  ...  0 

VW=  :  :  :  (25) 

0  0  ...  R{k)  _ 


D.  State  and  Covariance  Update 

The  state  update  equation  is  given  by 

x(k\k)  =  x(k\k  —  1) 

m{k)  ciprn(k) 

+W/V5,nv  (k)  E  E  P<P,nv{k)Vv,nv{k)  (26) 

(p= 0  nip  =  l 

and  the  covariance  associated  with  the  updated  state  is 

P{k\k)  = /30(k)P(k\k-l)  +  (l-/30(k))Pc{k\k)+P(k)  (27) 

where  the  covariance  of  the  state  updated  with  the  correct 
measurement  is 

Pc(k\k)  =  P(k\k-l)-W^nv{k)S^nv(k)W^nv(k)'  (28) 
and  the  spread  of  innovation  term,  P{k),  is  given  in  (29). 

E.  MD  Association  Probabilities 

There  will  be  m(k)  validated  measurements  at  time  k. 
Among  these  validated  measurements  one,  two  or  ip  number  of 
measurements  can  be  target  originated.  Multiple  detection  as¬ 
sociation  probabilities  are  evaluated  by  probabilistic  inference 
which  is  made  on 

•  number  of  measurements  in  the  validation  region,  m(k) 

•  number  of  target  originated  measurements,  ip 

•  location  of  measurements 
which  is  expressed  as 

P<p,nv(k)  =p{eVtnv(k)\Zk1m(k)1p,Zk~1)  (30) 
Applying  Bayes’  theorem 

Pv,nv{k)  =  ^p(Zk\6ip^v(k),m(k),p,  Z^1) 

xp(°^,nv(k)\m(k),p,Zk~1)  (31) 

The  first  term  in  (31)  refers  to  the  joint  density  of  the  pdf 
of  the  correct  measurement  is  given  in  (32)  where  Pq  is  the 
factor  that  accounts  for  restricting  the  normal  density  to  the 
validation  gate.  The  second  term  in  (31)  is  the  probability  of 
the  association  events  conditioned  only  on  rri(k)  and  p. 

lv,nv{k)  =  p(6(^nv(/c)|m(fc),v5,.Zfc_1) 

=  P{0<p,nv{k)\m(k),ip)  (33) 

where  the  probability  7v,n  evaluates  the  event  condi¬ 
tioned  on  the  total  number  of  validated  measurement  JA  =  to. 

Here  M  denotes  the  random  variable  and  m  its  realization  [2]. 

7 <p,nv(k)  =  p(Ov,nv\M  =  m(k),tp) 

=  p(^v,nv\^  =  m(k)  -  ip,M  =  m(k)) 
xpj'E'  =  m(k)  —  <p\M  =  m(k )) 

+p(8v,nv  I®  =  m(k),M  =  m(k)) 
xpj'f  =  m{k)\M  =  m(k ))  (34) 

where  T1  is  the  number  of  false  measurements.  For  ip  target 
originated  measurements,  T1  must  be  either  m{k)  —  p  or  m(k). 
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P(k)  =  Wv<nv{k) 


m(k )  c¥ 

E  E  ^nv{k)v^nv{k)vv^v(k)'  -  v(k)v(k)' 

ip= 0  n<p  =  l 


^W*)' 


(29) 


p(Zfc|^inv(/e),m(fc),^Zfc-1)a 


Pg 


xV(fc)  m(fc)+W(uViriv(fc);0,  S'(fc))  nv  =  1, ...,  ciprn(k) 


V(k)~ 


%(k) 


(32) 


=  0 


7  <p,nv(k)  —  \ 


where 


f  ^k)  x  pO  =  m(k)  -  ip\M  =  m(k)) 

Tl<p  —  1,  . . Ciprn  ( k ) 

p(\K  =  m(k)\M  =  m(k)) 


nv  =  0 


(35) 


p(\S?  =  m(fe)  —  ip\ M  =  m(k)) 

p(M  =  m(k)  |'f  =  m(k)  —  p)p{^  =  m(k)  —  p) 
p(M  =  m(k)) 

_  PDipPGp,(m[k )  -  p) 
p(M  =  m(k)) 


(36) 


and 


p(v]>  =  m(k)\M  =  m(k )) 

p(M  =  m(k)  |\P  =  m(k))p(\ D  =  m(fc)) 
p(.M  =  m(k)) 

=  (1  -  PDPG)p(m{k)) 
p(M  =  m(k)) 

Pov  is  the  probability  of  detecting  a  target  ip  times  per 
scan.  The  total  probability  of  detection  PG  will  become 
the  superposition  of  detection  probabilities  of  Pnv.  Also, 
Pd^Pg  is  the  probability  that  the  target  has  been  detected 
and  tp  measurements  originated  from  it  are  inside  the  gate 
and  (1  —  PnPG)  is  the  probability  that  the  measurements  in 
the  gate  are  false  alarms.  Thus 

m(k) 

p(M  =  m(k))  =  ^2  PDvPGp(m{k)  -  <p) 
v>— i 

+(1  -  PDPG)p(m(k))  (38) 


Substituting  (38)  in  (36)  and  (37),  the  result  in  (35) 

(  _ 1 _ PDvPgii(m(k)-<fi) _ 

m(k)  PDvPG^(m{k)-ip)+(l-PnPG)p{m{k)) 


1ip,nv  {k) 


< 


nv 


1, ...,  c(prn(k') 


_ {^-PpPg)k-(m(k)) _ 

PDvPgn{m(k)-ip)+(l-PDPG)ii(m(k)) 


nv  =  0 


(39) 


For  the  probability  mass  function  of  the  number  of  false 
measurements,  p(m(k)),  a  Poisson  or  diffused  prior  model 
can  be  used  in  the  volume  V(k)  (see  section  III-B). 

•  Poisson  model  (parametric  MD-PDA): 

M(m(fc))=e-^W(AT/(y)  (40) 

m{k)l 

where  A  is  spacial  density. 

•  Diffuse  prior  model  (non-parametric  MD-PDA): 

p(m{k))  =  fi(m{k)  —  p)  =  K  (41) 

where  K  is  a  constant. 

Finally,  by  substituting  (39)  in  (31),  the  association  proba¬ 
bilities  for  a  measurement  or  ip  measurements  set  can  be 
computed  with  parametric  (40)  or  non-parametric  (41)  false 
measurements  model  of  MD-PDA. 

IV.  Simulations 

A  surveillance  region  covering  an  area  of  1000  m  long  and 
1 500  m  wide  is  used  to  test  potential  advantages  of  the  multiple 
target  originated  measurement  approach.  Measurements  are 
generated  by  a  2D  radar  with  the  following  properties: 

•  Pd i  =  0.05  is  probability  of  detecting  a  target  once  per 
scan  of  the  measurement  data 

•  Pd2  =  0.9  is  probability  of  detecting  a  target  twice  per 
scan  of  the  measurement  data 

•  Pd  =  Pin  +  Pd2  =  0.95  which  is  total  probability  of 
detecting  a  target  in  a  scan  of  the  measurement  data  (i.e, 
Pd  used  for  PDA) 

•  Pfa  =  10-5/m2  with  Poisson  distribution 

A  single  target  starts  from  origin  and  moving  with  constant 
speed  of  15  m/s  parallel  to  the  x-axis  is  considered.  Target  ini¬ 
tialization  is  done  using  two  point  target  initialization  method. 
The  scan  interval  (  sampling  period)  is  1  s  and  it  consists  of  50 
scans.  For  the  MD-PDA  the  probabilities  of  detections  used 
are  Pm  and  Pd2  while  Pd  =  Pdi  +  Pd2  is  total  probability 
of  detecting  a  target  used  in  PDA. 

Figure  2  shows  the  Root  Mean  Square  Error  (RMSE) 
for  position  that  demonstrates  the  improved  performance  of 
multiple  detection  approach  over  the  classic  probability  data 
association.  As  PDA  tends  the  apportion  the  weight  among  the 
target  originated  measurements,  MD-PDA  assigns  the  weight 
to  measurement  set,  rather  than  a  single  measurement,  that  are 
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Figure  2.  Position  RMSE  evaluation  for  MD-PDA  vs.  PDA 


Time  (s) 


Figure  3.  Velocity  RMSE  evaluation  for  MD-PDA  vs.  PDA 

originated  form  a  target.  The  velocity  RMSE  evaluation  result 
is  shown  in  Figure  3. 

The  performance  evaluation  result  which  is  based  on  1000 
MonteCarlo  Runs  is  presented  in  Table  I.  With  respect  to 
Position  and  Velocity  RMSE  the  MD-PDA  performs  better 
than  PDA.  This  is  because  unlike  PDA,  the  MD-PDA  updates 
the  filter  with  the  set  of  measurements  that  are  originated  form 
a  target.  Due  to  more  association  events  evaluation  the  MD- 
PDA  takes  longer  time  than  PDA. 

Table  I 

Performance  Evaluation  (MD-PDA  vs.  PDA) 


Performance  Matrix 

MD-PDA 

PDA 

Position  RMSE 

1.62  m 

2.83  m 

Velocity  RMSE 

0.35  m/s 

0.74  m/s 

Average  Latency 

0.11  s 

0.08  s 

V.  Conclusions 


measurements.  When  multiple  detections  from  the  same  tar¬ 
get  fall  within  the  association  gate,  the  standard  PDA  filter 
returns  degraded  estimation  results  due  to  violation  of  one 
measurement  per  scan  assumption.  In  the  proposed  MD-PDA, 
combinatorial  association  events  are  formed  to  handle  the  pos¬ 
sibility  of  multiple  measurements  from  the  same  target.  Mod¬ 
ified  association  probabilities  are  calculated  with  the  explicit 
assumption  of  multiple  detections.  Experimental  results  show 
the  effectiveness  of  the  proposed  algorithm.  Similar  extensions 
to  handle  multiple  targets  using  the  JPDA  filter,  MHT  and 
MFA  tracker  are  under  development.  Also  further  work  has  to 
be  done  to  initialize  targets  with  multiple  detections. 
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In  this  paper  a  new  Multiple  Detection  Probabilistic  Data 
Association  (MD-PDA)  filter  was  proposed.  The  algorithm 
is  designed  for  tracking  a  target  while  receiving  multiple 
detections  form  the  same  target  within  the  same  scan  of 
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